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ABSTRACT 


We  consider  the  problem 


<p) 

< 

u  =  ( u  ) 
t  XX 

with 

x  e  R,  t  >  0 

u(x,0)  »  Uq(x) 

for 

X  €  R  , 

where  m  >  1  and  uQ  is  a  continuous,  nonnegative  function  that  vanishes 

outside  an  interval  (a,b)  and  such  that  (u®-1)"  <  -C  <  0  in  (a,b).  Using 

a  Trotter-Kato  formula  we  show  that  the  solution  conserves  the  concavity  in 
time:  for  every  t  >  0,  u(x,t)  vanishes  outside  an  interval  Jl(t)  “ 

(^(t)  J  and 


(um_1)  < 

'  XX 


C 

1  + 

m-1 


t 


in  Q(t).  Consequently  the  interfaces  x  -  CjU),  i  »  1,2,  are  concave 

curves.  These  results  also  give  precise  information  about  the  large  time 
behaviour  of  solutions  and  interfaces. 
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SIGNIFICANCE  AND  EXPLANATION 


The  flow  of  a  gas  through  a  porous  medium  is  governed  by  a  degenerate 
quasilinear  parabolic  equation*  It  is  known  that  the  nonnegative  solutions  to 


this  equation  possess  a  lower  bound  for  the  second  derivative  of  the  pressure 
in  the  spatial  variables,  This  bound  plays  an  important  role  in  the 


mathematical  treatment  and  is  related  to  the  entropy  of  the  flow.  Since  the 

Tob  X 

solutions  exhibit  interfaces  across  which  vA  jumps  positively,  no  upper 
bound  is  possible  globally  for  v,^.  Nevertheless  we-pr©ve  that  the  concavity 


of  vtj/t)  in  the  region  where  v  is  positive  is  preserved  in  time.  This  is 
in  itself  an  interesting  geometric  property  of  the  solution.  It  also  allows 
one  to  obtain  precise  information  about  the  asymptotic  behaviour  of  the  flow. 
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CONCAVITY  OF  SOLUTIONS  OF  THE  POROUS  MEDIUM  EQUATION 


Philippe  Benilan*'*^  and  Juan  Luis  Vazquez**'^  * '  ^ 

Introduction 

We  consider  the  Cauchy  problem 

f  u.  ■  (um)  in  8  «  R  <  (0,»l 


1 

f  u  -  (um) 

in 

(P)  l 

1  t  XX 

1 

l  u(x,0)  -  uQ(x) 

for 

where  m  >  1  and  uq  is  a  non-negative  continuous  function  vanishing  outside 
a  bounded  interval  (a,b)  with  Uq(x)  >  0  for  x  e  (a,b). 

It  is  well  known  (see  for  instance  [V2]  for  general  references  on  this 
problem)  that  there  exists  a  unique  generalized  solution  u(x,t)  of  (P)  that 
is  non-negative  bounded  and  continuous  on  Q.  Moreover  there  exist  two 
continuous  monotone  curves  x  *  (^(t),  x  *  £2^  called  the  interfaces  such 


that  £^(0)  -  a,  =  b  and 

a  d2f  { (x,t ) eQ;  u(x,t)  >  0>  -  {(x,t)cQf  ^(t)  <  x  <  ?2(t)} 

Also  u  €  C“(ft). 

In  the  description  of  the  flow  of  a  gas  through  a  porous  medium,  u 

represents  the  density  of  the  gas  and  v  ■  — ura_1  the  pressure.  We  shall 

m- 1 

retain  this  denomination  in  the  sequel.  Of  course  v  e  C"(fl)  and  it 
satisfies  on  Q  the  equation 

(0.1)  vt  -  (m-1)vvxx  +  v*  , 
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We  recall  also  the  basic  "convexity  inequalities* 


<°*2>  vxx  >  •  (bh*  1 ) t  111  D'(C) 

(0.3)  (-1)i(q  •*•  >  o  in  D'(*+> 

Now  we  assume  that  the  initial  pressure  Vq  ■  Uq”1  is  concave  in 
(a,b);  more  precisely  we  assume 


(0.3) 


(0.4) 


v0,xx  <  in  D' (a»b> 


with  some  constant  C  >  0.  The  purpose  of  this  paper  is  to  study  the 
concavity  of  the  pressure  v  on  Q  and  of  the  interfaces  on  R+,  a 

their  consequences.  The  main  result  is 


THEOREM  1.  Under  assumption  (0.4).  the  following  "concavity  inequalities" 


■xx  1+(m+1)Ct 


(0.6)  +  I.'(j5)ct  qi  <  °  1"  °M»+>  i2£.  i-1.2  • 

In  the  case  C  =  0,  the  concavity  (0.5)  has  already  been  proved  by 
Graveleau  and  Jamet  [GJ].  For  proving  (i)  in  the  general  case  C  >  0,  we 
will  follow  their  approach  based  on  splitting  (0.1)  into  the  two  equations 
(0.7)  vt  =  (m-1)wxx 

(0.8)  vt  -  v2  . 

We  then  recover  (0.1)  from  (0.7),  (0.8)  via  a  Ttotter-Kato  formula. 

In  the  case  C  >  0,  the  combination  of  (0.5)  and  (0.6)  with  (0.2)  and 
(0.3)  gives  strong  information  about  the  asymptotic  behaviour  of  the  pressure 
and  the  interfaces.  It  was  proved  in  (VI]  that  this  asymptotic  behaviour  is 
determined  in  first  approximation  by  the  two  invariants  of  the  motion,  namely 
the  total  mass 


(0.9) 


M  =  /  Uq ( x ) dx  =  J  u(x,t)dx 


and  the  center  of  mass 


(0.10)  x0  *  M_1  /  xu0(x)dx  -  M"1  /  xuQ(x)dx  . 

In  fact  there  is  a  selfsimilar  (Barenblatt)  solution  given  in  terms  of  the 
pressure  by 


(0.11) 

with 

(0.12) 

such  that  as  t  ♦  “, 
o(t_1)  uniformly  in 
(-1 )ir' (t)  +  o(t_1). 


v(x,t) 


Cr(t)2  -  (x-xq)2]+ 
2(o+1)t 


r(t)  *  c(rf““^t)  ,  c  ■  c(m)  >  0  , 

v-v  =  o(t_(m+1))  uniformly  in  x  c  R, 
n,  Ci(t)  =  xq  +  (-nMt)  +  o(1)  and 


vx  +  x/(m+1)t  * 
C[(t)  * 


As  a  consequence  of  Theorem  1  these  results  can  be  improved  as  follows: 
THEOREM  2.  Under  the  assumption  (0.4)  with  C  >  0  we  have  as  t  +  • 

(0.13)  vxx(xft)  a 


— ] -  +  0(-^~) 

(nri-1  )t  t2; 


uniformly  in  x  for  (x,t)  e  fl  and 

(0.14)  5j(t)  -  (-1  )ir"(t)  (1  +0 (•!))  . 


These  rates  of  convergence  are  optimal  as  it  can  be  easily  checked  on  the 
time-delayed  Barenblatt  solutions  v(x,t;x)  *  v(x,t+x),  where  v  is  given  by 
(0.11)  and  x  >  0.  From  (0.13),  (0.14)  sharp  estimates  follow  for  v,  vx,  5 


and  5 1 . 


The  main  problem  in  proving  concavity  results  lies  in  the  fact  that 
across  a  moving  interface  of  a  solution  to  problem  (P)  vx  has  a  jump 
discontinuity,  therefore  vxx  is  very  singular  and  positive  on  the 
interfaces.  Using  equations  (0.7),  (0.8)  allows  to  overcome  this  difficulty 
because  the  interfaces  in  (0.7)  do  not  move  while  in  (0.8)  we  are  reduced  to 
study  the  characteristic  lines  of  a  first-order  equation. 
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On  the  contrary  in  proving  lower  bounds  for  no  such  difficulty 

arises  since  we  can  always  approximate  our  solution  by  smooth  positive 
solutions  and  apply  the  maximum  principle  to  the  equation  satisfied  by  vH 
as  in  the  proof  of  (0.2)/  cf.  [AB].  Thus  if  v  is  a  solution  of  (0.1)  whose 
initial  datum  Vq  satisfies 

(0.15)  v0,xx  >  “C  in  D' (*) 

then 

(0.16)  vxx  >  -  jct  in  D' (Q) 

It  is  worth  noting  that  we  can  always  use  the  maximum  principle  method  of 

[AB]  when  dealing  with  problem  (P)  for  0  <  m  <  1  (the  so-called  fast- 

diffusion  case)  or  m  =  1  (the  heat  equation)  and  the  above  results  (0.5), 

(0.16)  are  true  if  we  define  v  *  -  —  u“^1_in^  if  m  <  1,  v  *  log(u)  if 

I  m 

m  =  1.  We  shall  leave  the  verification  of  these  facts  to  the  reader,  but  let 
us  recall  that  v  e  c"(Q)  in  both  cases  and  v  <  0  if  m  <  1. 

The  plan  of  the  paper  is  as  follows:  in  section  1  we  study  the  Cauchy 

problems  associated  to  (0.7)  and  (0.8)  for  vq  satisfying  (0.4).  We  prove 
the  statement  about  vxx  in  Theorem  1  via  the  Trotter-Kato  formula  in  section 
2.  Finally,  section  3  studies  the  interfaces  and  the  asymptotic  behavior. 


I 


a 


1 .  PRELIMINARIES 


In  this  section  we  study  the  Cauchy  problems  for  equations  vt  *  Xwxx 
2 

and  vt  *  vx<  To  make  their  application  in  Section  2  easier  it  is  convenient 
to  introduce  the  following  notations.  Let  N  be  the  set  of  continuous  non¬ 
negative  functions  f  :  R  +  R  satisfying  for  some  bounded  interval  1(f)  = 
(a,b) 

f ( x)  -  0  if  x  i  (a,b)  ,  f"  <  0  in  D'(a,b)  . 

It  is  clear  that  if  f  e  hi,  then  either  f  =  0  or  1(f)  =  {xeR;  f(x)  >  0}. 

W  is  a  closed  set  in  the  span  Cc(R)  of  continuous  real  functions  on  R 
with  compact  support  with  the  usual  topology.  To  measure  the  concavity  of  a 
function  f  e  N  we  define 


C(f)  -  sup{C  >  0  »  fxx  <  -C  in  D* (I (f ) ) } 
c(f)  =  sup{c  >  0  j  cfxx  >  -1  in  D'(R)} 


with  the  convention  C(f)  *  0  if  f  =  0.  Clearly  0  <C(f)<c(f)_1<«» 

and  the  maps  f*-*-C(f)/  f  »  elf)  are  upper  semi  continuous  functionals  on  N 
(endowed  with  the  topology  of  Cc(R)). 

Let  =  {fcN;  c(f)  >  0}.  Hie  following  results  can  be  proved  by 
elementary  calculus. 


,  1. 1. 

Let  f  £  N, 

with 

f  t  0  and  1(f) 

=  (e 

i) 

f  f  W2'°°(a,b) 

and 

C(f)  =  inf  ess  -  f" 
(a,b) 

<  sup  ess  -  f"  = 
(a,b) 

c  (f ) ' 

ii) 

f  e  W1,eo(R) 

and 

f  <  (b_a)2 

(f  ) 

,  2  R  f  K 

2  < 

f  <  c(f)  ' 

c(f ) 

iii) 

f'  e  BV(R)  , 

f’  (a 

+  )  >  0  >  f(b-) 

and 

-5- 
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f"  "  f"X(4,b)  +  f'(a+)5a  -  f*(b-)5b  in  D'(R). 


We  solve  first  the  Cauchy  problem 
( 1  •  4a )  vt  *  in  Q  ■  R  x  (Of08)  , 

(1.4b)  v(x,0)  ••  Vq (x)  for  x  e  R  , 

where  X  is  a  positive  constant  and  Vq  e  M.  We  obtain 

PROPOSITION  1.2.  Let  vQ  e  N  and  I(vq)  =  (a,b).  There  exists  a  unique 

function  v  e  C(Q)  that  solves  problem  (1.4)  in  the  sense  that 

(a)  v  >  0  _in  Q,  v  is  positive  precisely  on  the  set  G  *  (a,b)  x  (0,«) 
and  v  e  c“(£)) . 


(b)  vt  =  Xwxx  _in  fi. 


(c)  v(x,0)  =  Vq ( x)  for  every  x  e  R. 

Moreover,  if  v(t)  *  v(*,t)  then  for  every  t  >  0,  v(t)  e  and 
(1.5)  0  <  V(t)  <  Vq  ill  R  , 


(1.6) 


1 _  <  C<V 

c(vn)+Xt  *  V^xx  *  1+C(v„)Xt 


in  D'(a,b) 


If  c(v0)  >  0  then  for  every  (x,t)  e  Q 


(1.7) 


XtvQ(x) 

Vq ( x)  -  v( x, t )  <  c(^}  ■ 


Proof.  Let  us  assume  first  that  c(vq)  >  O.  Tb  avoid  the  degeneracy  of  the 


equation  (1.4a)  at  the  level  v  *=  O  we  begin  by  considering  the  problems 


vfc  =  >wxx  in  ft  =  {a  <  X  <  b,  t  >  0} 


v(x,0)  =  vQe(x)  for  a  <  x  <  b 


v(a,t)  =  v(b,t)  *  e(XtC{vQ)  +  1)"  for  t  >  0  , 


where  e  is  positive  and  {vQe}e  is  a  sequence  of  C*,  positive  functions 
in  [a,b]  such  that 

i}  ^v0e  )xx  <  c<v0>  and  <v0e}xx  >  Vc(v0)  , 
ii)  v0c  ♦  v0  uniformly  in  [a,bj  as  e  ♦  O  , 

iii)  vQe(x)  »  e  and  vQejXX  =  -C(v0)  at  x  «  a,b* 

Since  (P£)  is  a  quasilinear  parabolic  problem  and  we  have  the  a  priori 

estimate  e  <  v(x,t)  <  supx{vQe(x)}  for  the  solutions,  by  the  standard  theory 

(cf.  [LSU] )  there  exists  a  solution  vg  e  C°°(ft)  of  problem  (Pg). 

Consider  now  the  function  p(x,t}  =  vgxx(x,t).  It  satisfies  in  ft  the 

equation 

Pt  *  AvePxx  +  2XvexPx  +  *P2  ' 


with  initial  condition  p(x,0)  ■  Vgg  xx(x)  and  boundary  conditions 

C(V 

P(a,t)  -  p(b,t)  *  -  rr;  • 


r  XtC(vQ)+1 

(Note  that  p  *  vgt/(Xv).)  Since  for  every  k  >  0  the  function 

P(x,t;k)  ■  -(Xt+k)”^ 

is  an  explicit  solution  of  the  equation  satisfied  by  p,  it  follows  from  the 
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maximum  principle  that 


(1.8) 


Xt+c( Vg) 


<  p(x,t)  <  - 


C(v0) 


XtC(vQ)+1 


in  ft.  Since  p  =  v£fXX  and  v£(a,t)  =  v£(b,t)  <  e  it  follows  fromDl«nria 


1.1  and  (1.8)  that  the  following  estimates 


Xv 


>  v  >  -  '  - 

et  Xt+c(v. ) 


<  e  + 


( b-a ) ‘ 


2(Xt+c(vQ))  ' 


9  2,vn,» 

lv  I  <  - rrrr  and  |v  | 

1  ex1  c(v„ )+Xt  1  ex1 


(b-a) 


2(c(vQ)+2t) 


hold  in  (2.  Therefore  if  we  pass  to  the  limit  e  >0  there  is  a  subsequence 


{v  }  of  {v£}  that  converges  uniformly  in  [a,b]  x  [0,T]  for  every 


T  >  0  to  a  function  v  e  C (ft)  and  the  above  estimates  hold  for  v  so  that 


—  1  00 

v  e  W  '  ( (a,b)  x  (t,~))  for  every  t  >  0.  Because  of  (1.8)  in  a  neighborhood 


of  every  point  (x,t)  e  ft  the  sequence  v£(x,t)  is  uniformly  bounded  away 
from  0,  hence  by  standard  bootstrap  arguments  (cf.  [LSU]  or  [A])  we  conclude 


that  v  €  C  (ft)  and  vfc  =  Xv  vxx  is  satisfied  in  the  classical  sense.  It  is 


also  clear  that  v  e  c (ft)  and  both  the  initial  and  boundary  conditions  are 
satisfied. 

To  obtain  a  solution  v  in  Q  we  extend  v  by  v(x,t)  =0  if 
x  £  (a,b)  and  t  >  0. 

The  uniqueness  of  the  solution  is  a  consequence  of  the  maximum  principle 
(it  follows  in  particular  that  all  the  sequence  {v£}  converges  to  v) .  In 
fact  if  we  have  two  solutions  vfw  of  (1.4)  corresponding  to  initial  data 
vQ,  Wq  e  and  v0  >  w0  then  v  >  w  in  Q. 

In  case  c(vg)  **  0  we  can  yet  use  the  above  procedure  to  obtain  a 
solution  and  only  the  convergence  at  t  *  0  poses  a  problem  that  can  be 
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solved  by  approximating  the  initial  data  Vq  above  and  below  with  data  in 
and  using  the  above  comparison  results.  # 


REMARK.  Independently  of  the  authors  M.  Ughi  [U]  has  studied  the  equation 
^  =  uuxx  -  u(1-u)  in  R  x  (0,T) ,  for  initial  data  not  necessarily  concave. 
She  obtains  results  about  existence,  uniqueness  and  nonuniqueness  for  suitable 
classes  of  weak  solutions.  The  results  apply  to  equation  (1.4a)  with  obvious 
changes . 

Consider  now  the  initial-value  problem 


(1.9a) 


in  Q  , 


(1.9b) 


v(x,0)  m  Vq ( x)  for  x  e  R  . 


Equation  (1.9a)  is  one  of  the  simplest  and  best-known  examples  of  Hamilton- 
Jacobi  equation.  There  is  now  a  general  theory  of  viscosity  solutions  for 
equations  of  the  type  ut  +  H(ux)  *  0  where  H  is  a  continuous  function.  In 
particular  we  have  the  following  result  proved  by  Crandall  and  P.  L.  Lions. 
PROPOSITION  1.3.  [cf.  [CL],  Theorem  VT.2).  Let  v0  e  BUC  «/•  i  the  set  of 
bounded  uniformly  continuous  functions  in  R.  There  exists  a  unique 
v  e  C(Q),  which  is  a  viscosity  solution  of  ut  +  H(ux)  *  0  _in_  Q  with 
initial  data  v(x,0)  *  vQ(x). 

Moreover  the  maps  T(t)  ;  BUC(R)  ♦  BUC(R)  defined  for  t  >  0  _bj£ 

T(t) Vq  -  v(«,t)  form  a  strongly  continuous,  order-preserving  contraction 
semigroup  in  BUC(R). 

The  fact  that  in  problem  (1.9)  the  nonlinearity  H(p)  *  p2  is  convex 


considerably  simplies  the  construction  of  the  solutions  of  (1.9),  especially 


if  the  initial  data  Vq  r  W.  In  this  case  we  can  use  the  classical  method  of 
characteristics  (cf.  [Bt] ,  [L] )  as  follows.  Assume  for  simplicity  that 
Vq  e  C^Iajb)),  where  (a,b)  is  the  support  of  vQ.  Through  every  point 
(5,0)  with  a  <  5  <  b  we  construct  a  characteristic  line 
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(1.10) 


x(C,t)  -  5  -  2v$(£)t  . 


Along  this  line  the  derivative  w  =  vx  of  every  classical  solution  of  vt  = 
|vx|*  must  be  constant,  therefore  we  get 

(1.11)  V(x,t)  =  V(£,0)  +  /q  v(x(T),T)dT  »  V0(C)  -  |v^(5)|2t  • 

The  concavity  of  vQ  in  (a,b)  implies  that  these  characteristics  do 
not  cross  each  other.  In  that  way  a  C1  solution  of  (1.9a)  can  be 
constructed  in  the  region  where  v  >  0.  It  is  easy  to  see  that  this  region 
has  the  form 

(1*12)  Q  =  {(x,t)  e  Q  :  s^(t)  <  x  <  S2(tj} 

where  x  =  si(t)  and  x  =  S2(t)  are  Lipschitz-continuous  curves  (the 
interfaces),  s.j(0)  =  a,  S2(0)  “  b,  and  (-l)^s^(t)  is  nondecreasing  in 
time.  Along  the  interfaces  the  derivative  vx  is  discontinuous,  i.e.  a  chock 
propagates  according  to  the  well-known  Rankine-Hugoniot  conditions. 

The  following  results  about  problem  (1.9)  will  be  needed  in  the  proof  of 
Theorem  1 . 

PROPOSITION  1.4.  Let  v  be  the  viscosity  solution  of  (1.9)  with  v0  e  . 

Then  v  e  w1'°*(QT)  for  every  x  >  O,  v(t)  =  T(t)vQ  e  i  for  every  t  >  0 
and 

(1.13a)  c(v(t))  >  c ( vq )  +  2t  , 

c<v 

(1*13b)  C(v(t))  >  1+2tC( vQ )  ' 

Moreover  if  c(vg)  >  0  then 

2t#vnl,oo 

(1'14)  v<*>  -  v0  *  c(Vo)"  • 

Proof .  By  the  contracting  property  of  T(t)  we  may  assume  that  C(v0), 
c(vq)  >  0  and  vQ  e  C1([a,bJ).  To  prove  (1.13)  let  t  >  0  and  let  x1  and 
x2  be  two  points  s^(t)  <  x1  <  x2  <  s2(t).  There  are  well-defined  starting- 
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points  £.j  and  £2  (a,b)  for  the  characteristic  line  (1.10)  through 

(x1(t)  and  (x2,t)  respectively.  Setting  w1  -  vx(xvt),  w2  »  vx(x2,t)  we 
have 


x1  -  x2  *  Ci  -  C2  +  “  w1>  * 

But  vx  is  constant  along  the  characteristics,  hence  w2^*j 


since  vq  is  strictly  concave.  Therefore  if  we  write 


xrx2  5r52 


w2^i  W^, 


+  2t 


and  let  x^  -  x2  ♦  0  we  obtain  the  relationship 

1  _  1 


(1.15) 


+  2t 


vxx(x't)  v0(° 
valid  in  (2.  (1.13)  follows  immediately .  From  this  we  get 


(1.16) 


2IV-I 
2  c  0 

vt  "  vx  <  2t+c(v0)  * 


v^(x2)  -  vflfx-,)  >  0 


(1.14)  is  now  immediate. 


# 


TROTTER-KATO  FORMULA 


Let  us  now  consider  the  problem 
(2.1a)  vt  -  (m-1)wxx  +  vjj  in  Q 

(2.1b)  v(«,0)  -  v0  t  N  . 

To  solve  (3.1)  we  shall  use  a  Trotter-Kato  formula  based  on  the  results  of  the 
previous  section  (with  X  *  m-1).  It  is  interesting  in  that  respect  to  think 
of  the  solution  v(x,t)  of  problem  (1.4)  as  a  semigroup  8(t)  :  M  *  N,  t  >  0, 
defined  by  S(t)v0  *  v(t)  if  v(t)  *  v(«,t)»  v  being  the  solution  of  (1.4) 
with  initial  data  vQ.  In  the  same  way  the  solution  to  problem  (1.9)  defines 
another  semigroup  T(t)  :  N  *  SI,  t  >  0.  Both  are  order-preserving  continuous 
semigroups  in  SI.  The  concavity  estimates  (1.6),  (1.13)  can  be  reformulated 
as 

(2.2a)  c(S(t)v0)  >  c(v0)  +  (m-1)t  , 

(2.2b)  c(T(t)vg)  >  C(v0)  +  2t  , 

C(v  ) 

(2.2c)  C(S(t)v0)  >  1t(V,)C(,olt  ' 

and 

c<v0> 

( 2. 2d)  C(T(t)v0)  >  1+2C(v--)t  * 

Result  of  the  type  of  Chernoff  and  Trotter-Kato  formulas  for  pairs  of 
continuous  semigroups  are  known  when  both  semigroups  are  contractive  in  some 
Banach  space  X,  cf.  [BCP].  Unfortunately  in  our  case,  while  the  semigroup 
T  is  contractive  in  l“(B)  or  BUC(B),  the  semigroup  S  is  not.  In  fact 
the  only  contraction  properties  known  for  the  solutions  of  (1.4)  apply  to  the 
function  u  *  log(v),  which  satisfies  the  equation  ut  *  X(exp(u))xx,  cf. 
[BC]. 
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ttie  convergence  of  our  Trotter-Kato  formula  relies  on  the  compactness  due 
to  (2.2).  We  proceed  as  follows.  For  every  e  >  0  we  define  the  approximate 


solution  v£(t)  to  problem  (2.1)  by 


(2.3)  y  (t)  =  t-(n-1)e  v»  +  vn-1 

e  e  e  e  e 

if  (n-1)e  <  t  <  ne,  where  for  n  =  1,2,... 

(2.4)  v“  -  (S(E)T(e))n(v0)  , 

we  obtain  the  following  convergence  result 

PROPOSITION  2.  1.  As  e  0,  ve  converges  to  a  function  v  e  (Q) 
for  every  x  >  0  uniformly  on  compact  subsets  of  Q  =  R  x  [0,»), 

(2.5)  V(;,Stl(v2)xx-  ,m.2)v2 

in  D'(Q).  Moreover  if  v(t)  !  x  +  v(x,t)  then  for  every  t  >  0, 
and 


(2.6a) 


C(V 

C(v(t))  >  THfcnt  ' 


(2.6b) 


c(v(t))  >  c(v0)  +  (m+1)t  . 


Also 

(2.7) 


(2.8) 


2  .  2,V- 

x  c(v0)+(nrt-1)t  * 

(m-l)Jv  I  2*vn« 

0  oo  0  eo 

c(vQ)+(m+1)t  Vt  c(vQ)+(m+1)t 


Proof.  By  induction  we  have  vj?  c 


(2.9a) 

0  <  <  T(ne)v0  <  lvQlw  , 

(2.9b) 

_  C(v  ) 

e  1+(m+1 )neC(vQ)  ' 

(2.9c) 

c(v”)  >  c(vQ)  +  (nH-l)ne  , 

<2.9d) 

supp(v")  C  supp(T(ne) vq)  . 

n  w1  ,C8(qt  ) 

v  satisfies 


v(t)  e  1 


Observe  that  (2.9a)  follows  from  the  fact  that  S(e)v  <  v  and  T(e)  and 


S(e)  are  order-preserving  maps.  (2.9b)  and  (2.9c)  follow  from  (2.2)  and 
finally  (2.9d)  Is  a  consequence  of  (2.9a). 

Assume  now  that  c(vg)  >0  and  fix  tg  >  0.  We  have 
vg  e  W1#"(*  x  (0,tQ))  , 

0  <  ve  <  'V-  ' 

supp(vg(t))  C  supp(T(tQ+€ )vQ)  if  0  <  t  <  tQ  , 


v  <  2lvnl  /c(v„) 
e,x  0  •  0 


Moreover 


d+v  S(e)w  -w  T(e)v  -v 
_ e  _  e  e  _ e  e 

dt  e  e  ' 

where  v  (t)  *  v11  ^  if  (n-1)e  <  t  <  ne  and  w  *  T(e)v  .  Since 

C  £  C  t 

c(we(t))  >  C(v£(t))  +  2e  >  c ( Vg )  we  have  by  (1.7)  and  (1.14) 

(m-l)lvl  d+v  2lvl 

®  £  CO 

c(vQ)  <  dt  <  c(vQ) 

It  follows  in  particular  that  v£  is  relatively  compact  in  C(*  x  [0,tg]). 
If  ve  converges  uniformly  to  a  function  v  e  C(*  x  [0,tg]  )  along  a  sequence 
e  =  ek  •*>  0,  we  have  v  (  W1,-(R  x  [0,tg] ) 


0  <  v  <  T(t)v  ,  |v2|  < 

0  1  x1  c(v0) 

and 

(m-DIv  l.  2»v 

- . - Ji —  <  vt  <  r—  . 

c(vQ)  t  c(vQ) 

Furthermore  the  estimates  (2.6a)  and  (2.6b)  follow  from  (2.9b)  and  (2.9c)  and 
the  upper-semi-continuity  of  the  functionals  c  and  C. 

To  end  the  proof  of  the  proposition  we  have  yet  to  check  that  equation 
(2.5)  holds.  For  that  we  write 


-14 


(2.10) 


d  v 

ST  m  e  ^ols(T,"t't)),dT  +  7  * 

■  TT  /SKBttW.d))}^*,  ♦  J-  /S<(T(t);e(0)it)2dT  - 

-  Ell  f5((s,T)we(t))]t}2aT  a  I,  *I2  ♦  I,  . 

We  want  to  pass  to  the  limit  e  +  0.  Let  us  begin  for  instance  with  I2. 

Since  ve^  ♦  v  uniformly  in  R  x  [0,tg],  using  also  (1.14)  it  follows  that 

sup  *T(t)v  (t)  -  v( t )  I  =0  as  e.  ♦  0  . 

0<T<6  e  «  K 

0<t<tQ 

If  we  now  observe  that  c(T(x)v  (t) )  >  c(v0)  and  supp(T(x)v  ( t ) )  C 
supp  T(((t0+e)v0)  C  (-R, R]  for  some  R  >  0  the  convergence 
(2.11)  |  (T(T)v£(t))x|2  +  |v(t)x|2  in  CUO'tolrL^m)) 

will  hold  uniformly  in  x  e  (0,e)  as  a  consequence  of  the  following  compact¬ 
ness  result  applied  to  the  family  {f  (x)  •  T(x)v  (t)},  0  <  t  <  tn. 


LEMMA  3.2.  For  every  constant  c0,  R  >  0  the  set 

Nc0,R  *  ^  f  N  :  c(f)  >  c0,  supp(f)  C  [-R,R] } 


is  a  compact  subset  of 
constants  C  *  C(cq#R) 
(2.12) 

For  any  two  functions 


W1#P(R)  for  tuny  p  e  [1,*).  Moreover  there  exist 


>  0  and  6  c  (0,1)  such  that 

•frf2l  1  D  *  C,frf2,9n 

1  2  V1tP(R)  1  2  LP(R) 


f,f2  c  N 


C Qf  R‘ 


For  a  proof  of  this  result  see  [L,  Lemma  10.1].  It  is  based  on  proving 

that  N_  0  C  Ws,p(R)  for  some  s  >  1  and  p  >  1  and  then  using 
c0' K 

interpolation. 

Continuing  with  the  proof  of  Proposition  3.1  it  follows  from  (2.11)  that 
12  -  ~  fo{(T(x)vs(t))x}2  ♦  jv(t)x| 2  in  C((0,tQ]  ;  L1  (R) )  . 

In  the  same  way  we  prove  that 
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X1  S  e  /o{(s(T)we(t,)x}2dT  ■*  <v(t,x)2  in  c(l0'V  1  l1(m>) 


Finally  since 


Fg(t)  “7  /o(S(T)wg(t))2dT  +  v(t)2  in  C<*  x  [0,t0]) 

it  follows  that  I3  =  Fe(t)xx  ♦  (v(t)2)xx  in  D’<R  x  (0,t0)).  Therefore  we 
may  let  ♦  0  in  (2. 10)  to  obtain  (2.5). 


The  case  c{vQ)  =0  is  easily  dealt  with  by  comparison  since  for  every 


6  >  0 

i)  S(6)Vq  <  vQ  <  T(6)vfl  and  «T(6)vQ  -  S(6)v0l-  -*•  0  as  6  +  0, 

ii)  c(S(6)v0)  >  0,  c(T(6)v0)  > 

To  conclude  the  proof  of  the  first  part  of  Theorem  1  we  have  to  identify 
the  function  v  as  the  pressure  of  the  solution  to  problem  (0.1). 


PROPOSITION  2.3.  Let  v  (  C(Q)  D  W  '  (Q^ )  for  every  t  >  0  satisfy 
i)  for  every  t  >  0,  v(«,t)  e  N 


ii)  vt  *  (v2)xx  +  (2-m)|vx|2  in  D'(Q) 

iii)  v(.,0)  = 

where  un  is  nonnegative,  continuous  and  bounded.  Then  v  is  the  pressure 


associated  to  the  solution  of  (0.1). 


Proof .  Let  u  “  ((m-1)v/ra)  and  let  u  be  the  usual  solution  of  (0.1). 

By  the  preceding  results  if  0  is  the  subregion  of  Q  where  v  is  positive 
ft  has  the  form  {(x,t)  e  Q  :  s^(t)  <  x  <  s2(t)}  and  we  have  u  e  C (ft).  Also 
(by  standard  regularity  theory)  u  e  C"(ft)  and  ut  *  (um)xx  is  satisfies  in 
ft  in  the  classical  sense.  Since  u  <  u  on  the  parabolic  boundary  of  ft  by 
the  maximum  principle  we  get  u  <  u  in  ft. 

Moreover  since  (um)x  “  uvx  and  v(»,t)  e  N  it  is  easy  to  see  that 
<u»)x  +0  as  (x,t)  tends  to  the  lateral  boundary  of  ft,  hence  for  every 
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\  Vv', 


>  Ii' 


.  INTERFACES  AND  ASYMPTOTIC  BEHAVIOUR 


In  this  section  we  complete  the  proof  of  Theorem  1  by  studying  the 
concavity  of  the  interfaces  of  the  solutions  to  problem  (P)  and  derive  its 
consequences  for  the  asymptotic  behaviour  of  solutions  and  interfaces.  Theorem 


Proof  of  the  concavity  of  statement  (0.6)  of  Theorem  1. 

We  may  assume  that  v  is  a  solution  of  (0.4)  whose  initial  datum  vQ 
satisfies  (0.4)  with  C  >  0.  Also  we  may  consider  only  the  right-hand 
interface  and  drop  the  subindex,  ?(t)  =  ^(t).  Following  the  idea  of  the 
proof  of  the  converse  inequality  (0.3)  ([CF],  [VI])  we  compare  our  solution 
v  with  the  Barenblatt  solution  w  that  has  best  contact  with  it,  this  time 
from  above,  at  a  given  point  of  the  interface,  the  main  difference  lying  in 
the  fact  that  we  have  to  use  a  time-delayed  Barenblatt.  The  outline  of  the 
proof  is  as  follows;  let  tj  >  0.  Since  v(*,t1)  is  strictly  concave  in  its 
support,  cf.  (0.5),  vx(?(t1)  -  0,  t1)  is  necessarily  negative,  therefore  by 


the  interface  equation 
(3.  1) 


£'( t)  -  -v„(5(t)  -  0,  t) 


([A],  [K] )  we  have  £'(t).  We  consider  the  Barenblatt  solution 


(3.2) 


w(x,t) 


(r(t+x)2  -  (x-x1)2)+ 
~  2(m+1 ) (t+x ) 


where  r(t)  is  as  in  (0.12).  The  parameters  M,  t  >  0  and  x.,  e  R  shall  be 
adjusted  so  as  to  have 

(i)  v(C(ti),t.|)  =  w(C(t1),t1) 

(ii)  vx(?(t1)-0,t1)  =  wx(c(t1 )-0,t1 ) 

( iii)  v^U*,)  <  wxx(^,t1)  whenever  vt?,^)  >  0. 

By  (0.5)  and  (3.2)  the  last  inequality  is  implied  by 


1+(m+1 )Ct , 


(m+DIt^x)  ' 
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namely  t  >  1/((m+1)C).  The  best  choice  is 

t  *  1/{ (m+1 )C)  . 

We  can  now  choose  M  and  x<j  so  that  (i)  and  (ii)  are  satisfied.  It  follows 
from  (i)  -  (iii)  that  v(*,ti)  <  w(*,tj).  Since  both  v  and  w  are 
solutions  of  (0.7)  in  S  ■  R  x  (t^,*)  it  follows  that  v(x,t)  <  w(x,t)  in 
S.  Therefore  for  every  t  >  t^  we  have 

;(t)  <  Cw(t)  =  Xl  +  r(t+x )  . 

Since  at  t  *  t-|  we  have  5(t.|)  *  Cw(ti)  and  C'(t^)  -  C^(ti)  (*>» 

(ii))  and 

(3*3)  ?w(t)  “  "  (m+1 ) (t+T )  5w(t) 

we  conclude  as  in  [CF] ,  [VI]  that 

<3"*>  Ctt,1  *  ~  f(ti> 

in  the  sense  of  measures  in  (0,«),  i.e.  (0.6).  # 

(0.3)  and  (0.6)  together  imply  that  c"  is  a  locally  bounded  function. 
(0.6)  can  be  reformulated  as  saying  that  5' (t)/r' (t+r)  is  nonincreasing  as  a 
function  of  t. 

We  now  study  the  asymptotic  behaviour.  We  state  in  detail  the  results 
for  the  interface. 

PROPOSITION  3.1.  Let  v  be  a  solution  of  (0.4)  such  that  v0  satisfies 
(0.4)  with  c  >  0  and  let  t  -  1/((m+1)C).  Then  as  t  +  *  we  have 

(3.5)  r(t)  +  Xq  <  ?(t)  <  r(t+T )  +  xQ 

(3.6)  r'(t)  >  ;'(t)  >  r'(t+T)  . 

Moreover  ?(t)  -  r(t+x)  t  x0  and  5 * (t)/r' { t+r )  +  1. 

Proof.  Since  5 ' (t) /r' (t+T )  is  nonincreasing  and  £(t)/r(t)  ♦  1,  cf.  [VI], 
it  follows  that  (t)/r' (t+x)  +  1.  Therefore  £(t)  -  r(t+x)  is 
nondecreasing.  Since  r(t+x)  -  r(t)  ♦  0  and  ?(t)  *  r(t)  +  xg  as  t  +  * 
([VI])  we  have  ?(t)  -  r(t+x)  +  x0-  This  proves  the  right-hand 
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inequalities.  The  left-hand  inequalities  were  proved  in  [VI]. 

COROLLARY  3.2.  As  t  ♦  • 


(3.7) 

(3.8) 

(3.9) 


|?(t)  -  (r(t)  +  xQ)|  <  C^t 


m 

m+1 


|C'(t)  -  r’(t)|  <  C2xt 


|c"(t)  -  r"( t) I  <  C3Tt 


where  C^,  C2,  C3  depend  only  on  m. 

Proof .  (3.7)  and  (3.8)  follow  reap,  from  (3.5)  and  (3.6).  Then  the  estimate 

for  £”  follows  from  this,  (0.3)  and  (0.6). 

To  end  the  proof  of  Theorem  2  we  remark  that  (0.13)  comes  from  (0.2)  and 
(0.5).  By  integration  we  obtain  estimates  for  v  and  vx.  It  is  to  be 
remarked  that  the  estimates  for  v,  vx,  5  and  5’  are  valid  also  in  the  case 
of  symmetric  solutions  without  the  assumption  of  concavity,  cf.  [VI,  Theorem 
B] .  In  that  case  the  result  extends  to  several  space  dimensions.  Let  us 
finally  remark  another  consequence  the  above  results. 

COROLLARY  3.3.  The  interface  of  a  concave  solution  to  problem  (P)  has  zero 


waiting  time.  More  precisely  with  the  above  notations  we  have 
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